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ABSTRACT 


Title of Thesis: Electron Heating Processes in the Middle 

Ionosphere 

Tieh Chun Chang, Doctor of Philosophy, 1970 

Thesis directed by: R. T. Bettinger, Assistant Professor 

of Physics 

In this study, the equilibrium distribution functions for 

electrons in the ionospheric plasma have been calculated for 

the energy interval of 0 - 15 ev, utilizing data collected in 

a pulse probe experiment carried out by Dr. Philip T. Huang on 

board the sounding rocket NASA 18.12 which was launched from 

Wallops Island, Virginia into the normal daytime ionosphere. 

By using the densities and characteristic temperatures of 

superthermal and thermal electrons, the electron heat input 

3 3 

rates have been calculated to be between 3 x 10 ev/cm -sec 
4 3 

and 1.5 x 10 ev/cm -sec in the altitude range between ap- 
proximately 120 and 240 km. The correlation between the 
electron heat input rates and the heat loss rates in this 
altitude range has been studied. 

The equilibrium energy distributions of the secondary 
electrons in the middle ionosphere by solar ionizing radiation 
have been calculated for electrons with energies less than 
15 ev. A high energy electron distribution (E > 15 ev) con- 
tributes the retarded probe current which appears to be a 
linear function of the retarding probe potential. 
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CHAPTER I 


Introduction 

The ionosphere may be defined as the part of the 
earth's upper atmosphere where ions and electrons are pre- 
sent in quantities sufficient to affect the propagation of 
radio waves. It is considered to begin at an altitude of 
50 km and to extend outwards into the magnetosphere. The 
ionosphere is commonly divided into several layers, which 
are caused by the processes of photoionization, ionic re- 
actions and diffusion. The source of the ionizing mecha- 
nisms is a wide spectrum of solar X-ray and extreme ultra- 
violet (EUV) radiation. The spectrum consists of numerous 
emission lines generated in the chromosphere and corona of 
the sun, and varying amounts of contiruum radiations. 

The division of these layers is based on the observed 
altitude profile of the concentration of electrons as 
determined by radio wave propagation experiments. The 
various layers are identified by radio scientists through 
the use of the symbols D, E, FI, and F2. 

The D region is commonly located between altitudes of 
approximately 50 and 85 km. Its primary neutral constituents 
are N 2 and 0 2 , with traces of nitric oxide. N0 + , O* , and 
N 2 are the primary ions. In this region, the sources of 
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photoionization are believed to be Lyman alpha radiation 
and X-rays in the 1 to 8 & wavelength interval. The re- 
lative importance of X-rays in this range, however, is 
doubtful because the intensity of these rays is ordinarily 
quite low except during times of solar disturbances. An 
outstanding feature of D region is the existence of a 
large number of negative ions, caused by electrons attach- 
ing themselves to oxygen molecules to form the negative 
ion 0~ . The ratio of the negative ion density to electron 
density at 70 km is close to unity during the day. 

The altitude range between approximately 85 km and 150 
km has often been called the E region, and contains large 
amounts of and molecules. Above 120 km, the mole- 
cules have been dissociated to form atomic oxygen. EUV 

o 

and X-rays in the approximate range 8 - 140 A are usually 

considered as the major sources of ionization in this region. 

In addition, the photoelectrons produced by X-rays can 

produce secondary ionization in this layer. Monochromatic 

Lyman beta radiation here produces a layer of ions and 

electrons centered around 105 km with a noontime electron 

5 3 

density of approximately 10 /cm . In this region, electron 
concentration shows a strong diurnal effect with maximum 
occurring shortly after noon; the ionization, however, does 
not entirely disappear at night. Antonova and Ivanov- 
Kholodnyy (1961) have suggested that this night time 
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ionization is due to the influx of soft electrons from the 
magnetosphere . 

In the F region, the main neutral constituents are 0 

and N 2 . Traces of 0 2 , N, and NO are also common. The 

formation of this region is believed to be cuased by the 

extreme ultraviolet solar radiation in the wavelength 

range 150 to 900 A . The two maxima of electron density 

in the F region are known as the FI maximum, located around 

5 3 

200 km, with electron densities near 10 /cm , and the F2 
maximum, found between 250 and 400 km, with electron den- 
sities near 10^/ctcP . In this region, the peak in electron 
density distribution varies with the time of day, season, 
solar cycle, and latitude. Although the FI and F2 regions 
appear distinct, they are actually generated by the same 
ionizating radiation. Since the height of the FI peak 
(200 km) is different from the height of maximum produc- 
tion of electrons by solar radiation (150 km), it indicates 
that the formation of the F region is more complex than 
the explanations gives thusfar. 

The main production and loss processes in the F region 
can be expressed ass 


Production : 


A + hz^ 


■* A + + e 


( 1 . 1 ) 
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-A + + XY — » XY + + A 


( 1 . 2 ) 


Loss : 


XY + + e s- X' + Y' 


(1.3) 


where equation (1.1) represents the photoionization of an 
atom, A , by the ionizing radiation hz/ , equation (1.2) 
represents a charge exchange between the positive ion of 
A and a molecule XY, and equation (1.3) represents a dis- 
sociative recombination process. 

As the altitude increases, the effective recombination 
coefficient decreases, causing the electron density to 
increase beyond the peak of production. This may lead 
to the formation of a FI peak. Also at higher altitudes 
the chemical equilibrium of the F region tends to become 
a diffusive equilibrium. This process generally results 
in a F2 peak of electron density, with a maximum occurring 
during both day and night. 

Most of the structure of the ionosphere depends direct- 
ly upon the density and composition of the neutral atmo- 
sphere. Below 100 km, the main constituents of the at- 
mosphere are and 0 2 , with small amounts of Ar and C0 2 . 
Above 100 km, atmospheric mixing processes become less 
important than diffusion. Diffusive equilibrium may begin 
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at about 120 km, causing the general character of the at- 
mosphere to change. The density of each constituents 
decreases with height as given by the barometric law, can 
be expressed as; 


n ± (z, t) 


n. (z ) T < z o> exp(- 
T(z,t) 



dz/H i ) 


(1.4) 


where the scale height is defined as: 


H. 

l 


kT 

m ± g(z) 


(1.5) 


is the density of ith constituent with particle mass 
irn at height Z and at time t, T is the temperature, k the 
Boltzmann constant and g(Z) is the acceleration due to 
gravity at height Z. Each neutral gas tends to be dif- 
fusively separated and each density decreases with increased 
altitude according to its mass and temperature. At higher 
altitude, the gas temperature rapidly increases by a large 
factor, causing a slower decrease in the density of neutral 
atmosphere. 

In the bottomside inosphere, below the F2 maximum, 
photochemical processes strongly influence the composition 
of charged particles. The complex chain of events in the 
photoionization process during the day leads to the event- 
ual heating of the neutral atmosphere. Due to the relative- 
ly poor thermal contact between electrons and the more 
massive particles, a substantial temperature difference 
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will exist between the massive particles and electrons. 

Below 300 km, the assumption of a thermal equilibrium 
between ions and neutral particles is valid since at this 
height the ion gases are rapidly cooled by the neutral 
atmosphere. In this investigation, different aspects of 
the thermal structure of the middle ionosphere will be 
studied. Chapter II presents the results of rocket-borne 
pulse probes. In Chapter III the thermal electron heating 
and cooling processes are discussed. The photoelectron 
energy distribution function is discussed in Chapter IV. 
Summary of results will be presented in Chapter V. 
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CHAPTER II 

Pluse Probe Experiments 

In order to understand the thermal structure of the 
middle ionosphere, a number of experiments have been carried 
out in recent years. To date, the best opportunities to 
compare the theoretical understanding of the basic aspects 
of electron energy balance with the actual ionospheric 
environment have been provided by the results of rocket- 
borne Langmuir probes. The primary objective of these 
experiments was to measure the electron concentrations and 
energy distributions in the ionosphere. 

A. The Langmuir Probes 

The Langmuir probe is a collector, immersed in a 
plasma. The current is recorded as a function of the 
applied voltage and this information is analyzed in terms 
of concentrations and temperatures of the plasma consti- 
tuents. Mott-Smith and Langmuir (1926) were the first to 
deal extensively with the theory of probes with simple 
geometries including! (1) planar geometry, neglecting the 
edge effect; (2) cylindrical geometry, neglecting the end 
effect; and (3) spherical geometry. Different expressions 
for the probe current as a function of probe to plasma 
potential were derived for various particle velocity 


distribution. 
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The average speed of electrons is usually greater than 
that of ions. When a body or a probe is immersed in 
ionospheric plasma, it acquires a negative equilibrium 
potential. The probe repels negative ions and electrons, 
but attracts positive ions, and thus becomes surrounded 
by a positive "sheath" or region of positive space charges. 
The total positive charge is the sheath equals the negative 
charge on the probe. In their derivation, Mott-Smith and 
Langmuir avoided the electrostatic problem by assuming 
that the sheath surrounding the collector could be approxi- 
mated by one with a sharp outer edge, outside of which the 
plasma is not disturbed by the probe, and the potential is 
that of the undisturbed plasma. When the probe potential 
varies in relation to the plasma, the sheath thickness 
changes accordingly. Bettinger and Walker (1965) derived 
expressions for the sheath around the Langmuir probe for 
spherical and cylindrical geometry. 

In the use of a cylindrical probe, if the particles 
have a Maxwellian velocity distribution, the current can 
be expressed as the following equation (Mott-Smith and 
Langmuir, 1926) : 


(1 - erf ( 


R 2 

a 2 - R 


2 ^) + e 1 ^ erf ( 


a 2 ± 
a 2 - R 2 ' 


i 


27TR1 I 


a 

R 


'jr'?>o (2.1a) 
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= 2VR1 l o e^ 


< o 


(2.1b) 


where erf (x) 



dy 


complementary error function 


I 


O 



thermal current density 


f = 



= nondimensional probe potential 


R = the probe radius 


a = the sheath radius 
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In the accelerating region, the current collected by the 
probe is a function of the sheath radius, while in the 
retarding region the current is independent of the sheath 
radius. The dimensions of the sheath surrounding a 
cylindrical probe with a rather large ~jr are given in the 
following equation: (Bettinger and Walker, 1965) 


S 


( 2 / 3 )^ 

In (s/p ) 


\ 


h 



( 2 . 2 ) 


a _ R 

where S = — r — = nondimensional sheath thickness 


= = nondimensional probe radius 


h = = Debye length 

Y 47Tne^ 


When a probe is immersed in a plasma, the ratio of 
the sheath radius to the probe radius can be derived by 

0 

using equation (2.2). For a plasma with a density of 10 
electrons/crn and a temperature of 2000° K, the ratio of a/R 
has been calculated for a different applied potential . 
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The results are presented in Table 2.1 and plotted in 
Figure 2.1. 

For a Maxwellian energy distribution, the electron 
temperature may be measured by using the i-V characteristic 
in the retarding region. By taking the derivative of the 
logarithmic current with respect to the applied voltage, 
equation (2.1b) can be rewritten as: 


d log i _ e 

d V kT (2.3) 

hr 


A plot of log i as a function of the potential V 

P 

will produce a straight line whose slope is inversely 
proportional to the temperature T. Once the electron 
temperature is known, the electron density can be obtained 
by using the equation (2.1b) together with a knowledge of 
the vehicle potential. 

For a probe with a negative potential relative to 

the plasma, ambient electrons with energies E ^ eV and 

P 

with proper directions of incidence may reach the probe 
surface and contribute to the current. The electron 
velocity distribution for an isotropic plasma can be found 
by taking the second derivative of this current with respect 



TABEL 2.1 


Ratio of Sheath Radius to Probe Radius as a Function of Applied Voltage 


Applied Retarding Ratio of Sheath Radius to Probe Radius (a/R) 


Voltage (v ) 
P 

R = li 

R = 2h 

R = 3h 

i 

5.7851 

3.9008 

3.3100 

2 

8.1754 

5.2012 

4.1375 

3 

10.1663 

6.2850 

4.8995 

4 

11.9437 

7.2526 

5.5760 

5 

13.5789 

8.1413 

6.1962 

6 

15.1054 

8.9645 

6.7763 

7 

16.5496 

9.7430 

7.3230 

8 

17.9237 

10.4906 

7.8449 

9 

19.2510 

11.2049 

8.3425 

10 

20.5270 

11.8948 

8.8227 

11 

21.7642 

12.5558 

9.2862 

12 

22.9577 

13.2008 

9.7344 

13 

24.1238 

13.8286 

10.1753 

14 

25.2618 

14.4418 

10.6020 

15 

26.3756 

15.0352 

11.0187 


R — probe radius 
h - Debye length » 



Assume T = 2000 K and n ~ 10*Vcm^ 



RATIO OF a/R 
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to the retarding probe potential (Druyvesteyn, 1930). 

For a cylindrical probe in an isotropic velocity distri- 
bution? neglecting the end effects, the Druyvesteyn 
relationship can be written as (Huang, 1969); 


dV ? 


P 


Ae~ 


2(2meV ) 2 


F(eV p ) 


(2.4) 


where A is the total surface area of the probe, and F(E) 
the distribution function satisfying the normalization 
condition: 


r 00 

F(E) dE = n (2.5) 

J 

o 


B. Gridded Probes 

One type of modified Langmuir probes are known as 
gridded probes, commonly referred to as "retarding potential 
analyzers." A number of gridded probes have been carried 
on various flight vehicles to determine the density and 
energy distribution of the ionospheric plasma (Hinteregger , 
I960? Sagalyn, et al. , 1963? Hanson, et al., 1964? Bettinger, 
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1964; Knudsen and Sharp, 1967; Harris, et al., 1967; Shea, 
et al., 1968; Moss and Hyman, 1968; and Huang, 1969). 

Of several types of gridded probes, the main one used 
in experiments referred to in this paper is the pulse probe. 
This was first developed and tested in flight by Bettinger 
in 1964. It consists of an outer cage and two closely 
spaced grids surrounding the collector. In operation, 
the entire probe is immersed in a plasma and a linear 
potential sweep is applied to it with respect to the vehicle. 
The pulse probe is operated in two modes. They are: (1) 

DC mode, and (2) pulse mode. The chief advantage of using 
the pulse techniques is that the electron concentration 
measurements do not depend on the velocity distribution of 
the electrons. The pulse current I shows a maximum at 

p 

the plasma potential and yields the ambient density without 
further recourse to a knowledge of the vehicle to plasma 
potential. 

When a probe is operated in the DC mode, ambient 
electrons with energies E > eV^ and the proper directions 
of incidence reach the surface of the collector. By letting 
Vp be the radial and v the tangential component of the 
velocity of an electron, with Vp defined to be positive 
when it is directed toward the origin, the total number of 
electrons which come to the collector in unit time may be 


written as: 
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i = 2ja 1 ■ e 


j 

—03 


V p f(v, v p ) dvdv^ 


( 2 . 6 ) 


where f(v,v^) is the ambient electron velocity distribution 
function normalized to the electron density. 


<33 

f r 


f (v, v^, ) dvdv^ = n 


(2.7) 


In the equation (2.6), a is chosen as the sheath radius, 
such that outside the sheath the distribution is not distrubed, 
and c< specifies the fraction of current incident on the 
sheath edge that actually reaches the collector, consistent 
with conservation of angular momentum and energy. 

The laws of conservation of energy and angular momentum 
are reconciled in the following equation: 


%m(v (R) + v^ (R) ) + 


eV = %m (v (a) 

kr 


+ V 


r 


(a) ) 


( 2 . 8 ) 
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v(R)R = v(a)a 


(2.9) 


where m is the mass of electrons and denotes the retarding 
probe potential. By using equation ( 2 . 9 ), equation (2.8) 
may be rewritten as: 


o o o -,2 2eV 

(R) = Vp (a) + v 2 (a) (1 - ~) £ 

r r r z m 


( 2 . 10 ) 


For a particle to reach the collector, the radial velocity 
should be at least greater or equal to zero at the probe 
surface, i.e., v^(R) a 0, or 


( a ) - v 2 (a) 



1 ) 


+ 


2eV 

£ 

m 



( 2 . 11 ) 


For a Maxwellian velocity distribution, the total current 
taken by the collector for a retarding probe is 
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i 


2vR 1 n e 


kT 

2^m 


e 


eV 


“ E 

kT 


( 2 . 12 ) 


When the probe is operated in the pulse mode, electrons 
can be accelerated toward the collector within the duration 
of the pulse. The electrons outside the cage near the probe 
are repelled away from the vicinity of the probe. The 
effect of the pulse on the motion of ions usually is 
negligible because of their relatively large inertia. 

The density of the electron n inside the cage can be 
expressed as a linear function of the pulse current I 


n 


fi £ f v 


eff 


( 2 . 13 ) 


where v ^ is the effective probe volume, f is the pulse 
repetition rate, and P is the combined transparency of the 
grids surrounding the collector. The density n inside the 
effective volume of the probe is equal to the density 
immediately outside the cage surface multiplied by the 
transmission coefficient of the cage which is assumed to be 0.8. 
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The pulse probe experiments developed by Bettinger in 
1964 were continued by Huang in 1969. In his investigation, 
Huang used a pulse probe consisting of a central collector 
surrounded by two closely spaced grids and an outer cage. 

The collector was plated with rhodium and had a diameter of 
approximately 0.8 in. The inner screen grids were made of 
stainless wires, and were located approximately 0.55 and 
0.65 in., respectively, from the axis of the probe. The 
outer cage was supported by three insulated circular rings 
providing the probe with overall dimensions of 2.9 in. in 
height and 2.8 in. in diameter. The pulse probe was mounted 
on top of the payload section (Figure 2.2). 

In the DC mode of operation, the two outermost 
elements of the probe were maintained at the same potential. 

A potential sweep was applied to the probe, varying from 
approximately plus 3 volts to minus 13 volts; the sweep 
rate was 54 volts/sec. The vehicle to plasma potential 
was approximately 2 volts throughout the flight and was 
obtained by the following method: first, extend the retarded 

thermal current which appeared as approximately a straight 
line on the semilog plot, towards more positive potential. 
Then, computing for the drift current I c by: 


I = 


< kT 
27rm 


D 


2ttri ne 


(2.14) 




Fig. 2.2 N/T 18.12 PAYLOAD SECTION 
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where n is the ionosonde electron concentration. Next, 
locate the point on the extended current curve corresponding 
to I Q . The applied potential corresponding to. this point 
was then equal to the negative vehicle potential. The values 
for the vehicle potentials obtained in this manner as a 
function of altitude are plotted in Figure 2.3. Under this 
arrangement of potentials, ambient electrons with sufficient 
energy to surmount the potential barrier and with proper 
direction of incidence were collected. 

In the pulse mode of operation, a negative square 
wave of 20 volts in amplitude was applied to the outer cage 
relative to the inner grid. A fast rising negative pulse 
was then applied to the outer cage. The pulse duration 
was approximately 0.2 / seconds, with a frequency of 100 
KHz. When the pulse was applied to the cage, all electrons 

inside the probe were accelerated toward the collector. 

-9 

It took only 5 x 10 seconds for an electron, initially 
at rest, to travel to the collector when plus 20 volts 
were applied to the outer cage. 

A pulse duration of 0.2/ seconds provided sufficient 
time for most of the electrons to reach the collector. 

Since electrons have high mobility, the plasma quickly 
reverts to a state of equilibrium, to which a new pulse 
may be applied. The effect of the pulse on the motion 
of ions, for the most part, can be considered negligible 



ALTITUDE (km) 
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because of their relatively large inertia. 

If the electron energy has a pure Maxwellian distri- 
bution, a plot consisting of log i versus applied retarding 
voltage will produce a straight line whose slope is pro- 
portional to the temperature. If the result is not a straight 
line, then the energy distribution is non -Maxwellian and 
can be distinguished from the results. 

C. Collection of Data 

In the DC mode of operation, the pulse probe current 
in a retarding probe potential has essentially the same 
current-voltage characteristic as the Langmuir probe except 
that there is no positive ion contribution to the total 
collector current. In this experiment, the collector current 
was monitored by a fixed range electrometer. The sensiti- 
vity of the electrometer was controlled by the value of the 
feedback resistor R, located between the input and the output 
terminals of the operational amplifier. The value at the 
output of operational amplifier corresponding to electrometer 
saturation was approximately 5.7 volts. Thus, the sensitivity 
range of the electrometer for a particular feedback resistence 
of value R (~10 5 , 10 8 , 1C)" 7 , 10 8 , 10 8 /1 ), was from zero to 
(5.7/R) amperes. The collected current corresponding to 
the maximum applied retarding potential was due to the 
accumulation of ambient electrons with energies greater 
than approximately 15.5 ev. By subtracting this quantity 
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from the ambient electrons with energies less than 15 ev. 

Since two distinct linear regions appeared when the 
current was plotted on a logarithmic scale, the current 
can be considered as arising from two distinct distributions, 
a thermal and a secondary distribution. Section AB shown 
in Figure 2.4 may be identified as the retarded thermal 
current and section BC as the retarded secondary current. 

The retarded thermal current produces an approximately 
straight line, and the majority of the electrons that 
contribute to this current may be characterized as a 
Maxwellian distribution. 

The results of the secondary probe current in this 
experiment approximated an exponential function (usually 
associated with a Maxwellian distribution of energies) 
of the retarding probe voltage over a rather extensive 
portion of the curve. In the experiment the exponential 
slopes resulting from the secondary probe currents were 
about 10 times smaller than those of the exponential 
slopes of the thermal current section. Further expansion 
of the probe current is expressed in Figure 2.5 which 
shows a linear section between approximately 6 to 13.5 


volts on the i-V plot. Since 
range roughly from 8 to 15 ev. 


~ 0 in the energy 
dv 

it may be concluded that 


there exist few electrons in this range, according to the 
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Druyvesteyn relationship. This linear current is attributed 

to ambient electrons with energies greater than ^ 15 ev. 

It is unlikely that this current is due to photoelectric 

or secondary emission electrons from the collector. The 

large positive potential bias on the collector 26 ev) 

relative to the inner grid would suppress electron emissions 

from the collector significantly. Displacement currents 

associated with the probe sweep were less than 2x1 0 ' L0 

amperes in this experiment. The displacement current 

was proportional to » the time rate of change for the 

dV 

sweep potential. In this experiment, ^ was approximately 
a constant. Therefore, the overall effect of the displace- 
ment current is approximately a constant terra. 

In the pulse mode of operation, all electrons inside 
the probe are accelerated toward the collector when the 
pulse is on. The density of the ambient electrons in the 
ionosphere may be derived from the pulse current by 
use of equation (2„13). Unfortunately, due to failures 
of the probe erection mechanisms in the experiment, only 
a fraction of the pulse probe had proper exposure and an 
adequate distance from the vehicle. As a consequence, the 
pulse probe was limited to retarding electron energy dis- 
tribution measurements. 
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CHAPTER III 

Thermal Structure of the Middle Ionosphere 
A. Introduction 

Measurements made by rockets and radars confirm a 
deviation in the thermal equilibrium near the 120 km level 
and show the rapid increase of electron temperatures in 
comparison with the temperature of ions and neutral 
constituents at an altitude equivalent to that of the F2 
peak. Several investigations have been carried out recently 
(Hanson and Johnson, 1961? Hanson, 1963? and Dalgarno, 
et al. 1963) to study thermal equilibrium within the 
middle ionosphere. 

Because of the small electron-to-ion mass ratio, 
the amount of heat imparted to the ion gas by superthermal 
electrons is small compared to that imparted to the thermal 
electron gas. The electron-to-neutral-particle mass ratio 
is also very small? therefore, in the elastic collisions 
the energy transfer between electrons is much more efficient 
than the energy transfer between electrons and ions or 
electrons and neutral particles. This makes the temperature 
of the electron gas, in general, higher than that of the 
ion gas or of the neutral constituents in the atmosphere. 

The superthermal electrons measured in this experiment 
may have been composed mainly of photoelectrons and secondary 
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electrons produced by high energy primary photoelectrons. 

The contribution which comes from a secondary distribution 
(apart from the thermal electron distribution) consists 
of electrons with energies ranging from approximately 
1 ev to 8 ev. For these electrons. Coulomb collisions 
with thermal electrons and inelastic collisions with neutral 
constituents are the dominant mechanisms of energy loss. 

The loss of superthermal electrons to ions is small compared 
with other losses. 

Hanson (1963) has indicated that up to an altitude of 
about 300km the photoelectron mean free path is shorter 
than the scale height of both the neutral and ionized 
components in the middle ionosphere. Hence, the rate of 
cooling of the superthermal electrons is assumed locally 
to be equal to the rate of heating at all altitudes (120 
km to 240 km). The contribution of the thermal conduction 
flux to the heat balance in the altitude range from 120 km 
to 240 km has not been considered in this investigation. 

The heat input rate from superthermal electrons to 
thermal electrons, has been calculated in Section B of this 
chapter. The cooling rates of thermal electrons produced 
by colliding with neutral constituents and encountering 
ions are discussed in Section C, and the results are pre- 
sented in Section D. 

B. Heat Input from Superthermal Electrons 
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In the NASA 18.12 sounding rocket flight, launched in 

1967, Huang measured the temperature of superthermal electrons 

(electrons with energies greater than ^1 ev) . The thermal 

electron temperatures were also measured in this experiment. 

The thermal electron densities were taken from the ionosonde 

electron concentration supplied by J. W. Wright of the 

National Bureau of Standards and are shown in Figure 3.1. 

Table 3.1 lists the density and temperature of both thermal 

electrons and superthermal electrons. At altitudes between 

150 and 240 km, the measured concentration of superthermal 

electrons is approximately constant, — 150 - 200 electrons/ 

3 

cm , and the equivalent temperature of the secondary dis- 
tribution has a range of 10,000 K to 30, 000° K. 

Since the temperature of the superthermal electrons 
is much higher than that of the thermal electrons, it 
will cause them to heat the ambient electron gas by elastic 
Coulomb collisions. The energy loss of a charged particle 
in a plasma by Coulomb collisions is expressed in a formula 
derived by Butler and Buckingham (1962). They derived the 
following equation by introducing the injection of a fast 
charged particle into a thermal plasma: 

r- 4 

dt = " Vwf - p ( u / w t ) ln A n e n s ev / cm3 ~ sec (3.1) 
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TABEL 3.1: CONCENTRATIONS AND TEMPERATURES OF 

THERMAL ELECTRONS /AND SUPERTHERMAL ELECTRONS 


Ascent 


Descent 


Altitude 

n 


T 


n 


T 


(km) 

5 3 

(lOVcin) 

(°K) 

2 1 
(10 /cni) 

(°K) 

120 

1.78 


5.8 

7,980 

126 

1.78 

1,422 

8.1 

8,117 

138 

2.11 

1,199 

3.0 

12,616 

149 

2.57 

1,390 

2.4 

12,508 

164 

3.11 

1,435 

1.5 

24,243 

174 

3.58 

1,594 

1.4 

23,308 

187 

4.24 

1,769 

1.4 

26,170 

194 

4.75 

1,978 

1.6 

26,767 

202 

5.17 


1.4 

31.328 

211 

6.01 

2,239 

1.4 

30,553 

217 

6.46 

2,305 

1.6 

27,721 

224 

7.31 

2,222 

1.6 

25,607 

228 

7.70 

2,165 

1.5 

30,820 

231 

8.13 


1.8 

27,907 

235 

8.40 

2,308 

1.9 

28,039 

237 

8.80 

2,338 

2.1 

29,261 

239 

9.00 

2,382 

2.4 

30,308 

237 

8.80 

2,291 

2.1 

28,241 

231 

8.13 

2,136 

1.6 

33,074 

227 

7.60 

2,065 

1.4 

32,601 

223 

7.20 

2,159 

1.5 

32.143 

216 

6.34 

2,029 

1.5 

29,642 

210 

5.85 

2,123 

1.4 

27,170 

204 

5.32 

1,924 

1.5 

26,426 

193 

4.70 

1,979 

1.5 

31,318 

184 

4.10 

1,804 

1.5 

27,334 

171 

3.50 

1,757 

1.5 

22,484 

162 

3.08 

1,461 

1.3 

23,465 

151 

2.73 

1,468 

1.8 

19,095 

134 

1.90 


3.6 

10,050 

123 

1.78 


11.5 

10,901 


In A 


15.42 
15.08 
15.20 
15.16 
15.24 

15.32 

15.43 


15.49 

15.50 
15.39 
15.32 


15.37 

15.37 

15.39 

15.34 
15.27 
15.26 

15.35 

15.32 

15.43 

15.32 

15.43 

15.36 

15.40 
15.19 
15.26 


Notes n g is the ionosonde electron concentration recorded during the time of 

Nike-Tomahawk 18.12 Launch (March 30, 1967), provided by J. W. Wright of 
The National Bureau of Standards. 
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where F(x) 


-x 


1 

X 


exp(-x 2 )dx - 2 exp(-x 2 ) 


and A is the ratio of Debye length to the distance of 
closest approach to the target in a head on collision 
(impact parameter equals to zero) 


3(kT e ) 3/2 

(47rn e )^e 3 


W t = thermal speed of the ambient electrons 
m = mass of electron, u = speed of test particle 

Another expression for the energy exchange between 
two Maxwellian gases at different temperatures was derived 
by Desloge in 1962. It can be expressed as: 


dU 

dt 


T 

8 n n In A (1 - =—•) 
e s e 

(27rmkT )** (1 + ^s) 3//2 
e T 

e 


ev/crtf 


sec 


where = Temperature of thermal electrons 

T g = Temperature of superthermal electrons 


(3.2) 
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The author calculated the energy loss rate for super- 
thermal electrons to thermal electrons through Coulomb 
collisions by using the above two equations. In these 
calculations, the approximate value of In A was 15, and 
varied slowly with electron density and temperature, 
therefore, it was assumed to be constant. Values of 
calculated In A as a function of altitude were tabulated 
in Table 3.1. These calculations show that Desloge's 
formula (3.2) constantly gives a value 30 to 40 percents 
higher than that derived from Butler and Buckingham's 
equation (3.1). The author's numerical results are listed 
in Table 3.2 and are depicted graphically in Figures 3.2 
and 3.3. 

C. Cooling Processes of Electron Gas 
1. Cooling by Ion Gas 

In the E and F regions the ions produced by photo- 
ionization are O , and C^. The relative production 

rates of these ions vary with the altitude in a manner 
determined by the composition of the neutral atmosphere 
and the atmospheric absorption of specific spectral ranges 
of solar radiation. NO + can be produced by either a 
charge exchange or ion-atom interchange reactions. In the 
daytime ionosphere, NO + and O* are the dominant ions 
below 165 km. Above this altitude 0 + starts to be dominant. 
At night, the transition level between atomic and molecular 



Ascent 
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TABLE 3.2s ELECTRON HEATING RATES IN THE MIDDLE IONOSPHERE 


Altitude 

(km) 


Butler and Buckingham Formula 
(10 ev/cm - sec) 


Desloge Formula 
(10 ev/cm - sec) 


126 

16.1 

19.4 

138 

5.8 

3.4 

149 

5.6 

7.9 

164 

3.1 

4.9 

174 

3.4 

5.2 

187 

3.3 

5.9 

194 

4.8 

7.3 

211 

4.9 

7.6 

217 

6.4 

9.5 

224 

7.5 

11.1 

228 

6.7 

10.4 

235 

9.3 

14.7 

237 

10.1 

16.7 

239 

12.7 

19.3 

151 

3.6 

5.6 

162 

2.7 

4.2 

171 

3.6 

5.4 

184 

3.8 

6.0 

193 

4.1 

6.4 

204 

5.0 

7.7 

210 

4.9 

7.4 

216 

5.7 

8.8 

223 

6.2 

9.6 

227 

6.0 

9.5 

231 

7.3 

11.5 

235 

9.6 

14.7 

237 

11.3 

17.0 

239 

12.3 

18.8 



Fig. 3.2 

ELECTRON HEATING RATE 
BY BUTLER a BUCKINGHAM METHOD 
. UPLEG 
4 DOWNLEG 
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ion dominance moves up to about 220 km. Other ions such 

* 4 “ *)■ 

as N , H^O are minor constituents. 

Hanson and Johnson (1961) have pointed out that 
Coulomb collisions with positive ions play an important role 
in the cooling of the electron gas. The energy transfer 
rate between two charged gases having Maxwellian velocity 
distributions with different temperatures can be written 
as s 


dU 

dt 


(ZjZ^e 2 ) 2 

’ 4 ^ n l n 2 m^m, - k < T l 


T ) 

1 2 


In A 


kT kT„ 3/2 

( — ~ + — -) 
m l m 2 


(3.3) 


Equation (3.3) is a generalized form of equation 
(3.2). For single charged ions colliding with electrons, 
equation (3; 3) is reduced to 


dU 


e 


dt 


i 


% 

m 2 

-4 /2TTn^ n. — k(T 
e i m. ' e 


T i) 


e 4 In A 

(kT ) 3</2 
e 


( 3 . 4 ) 


Numerically, taking In A Ct 15, it becomes 
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au . , (T e - V 

— — — — = -7.7 x 10 n n. r- 7 - 

dt e x T 3/2 

i e 


ev/cm - sec 


(3.5) 


where A. is the ion atomic mass in amu. 

1 

By comparing the coefficient with several 
different experimental measurements (Banks, 1966), the 
energy change rate may vary by + 10 % depending upon different 
ionospheric conditions. In the author's calculations, 
it is assumed that the positive ions are mainly 0 + . Under 
this assumption a small deviation will be introduced in 
the altitudes below 200 km. It is because that below 
200 km the energy transfer rate is dominated by collisions 
with neutral particles, while above 200 km the cooling by 
ion gas becomes the dominating heat loss mechanism for 
electrons. 

2. Cooling by the Neutral Atmosphere 

The three primary constituents of the neutral 

atmosphere ( 0 , N 2 ) are all effective in cooling the 

electron gas. The energy transfer rate ~ at which 

at 

unit volumes of electron gas convert their thermal energy 
to heat up these neutral constituents can be written as 
the sum of five terms: 
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dU 

dt 


(— ) 
l dt J o 


( dU) 

v dt ; N n 


(— ) 
v dt ; O, 


/ dU , 

dt) o A * 


( — ) 
'‘dt'Q 


M 


(3.6) 


where the first three terms on the right hand side of 
equation ( 3 . 6 ) are the elastic collision losses, (^ 7 ) 
is the atomic oxygen fine structure loss, and (4^) 


dt 0 


A 


dt O 


M 


is the atomic oxygen metastable state loss. Each of these 
terms is in general a function of electron temperature 
T , the neutral gas temperature T^, and the product of the 
electron density n g and the number density of each neutral 
constituent (^n(O), nC^), n (°20* The rotational loss of 
^ and C>2 is not very important above 120 km. 

The model atmosphere selected for the author's 
calculations was based on the neutral density profiles 
given by Hedin and Nier (1966). The data were measured 
in a rocket-borne mass spectrometer experiment. The neutral 
densities were then normalized to the neutral density adopted 
by Dalgarno, et al. (1963) at 120 km. The normalization 
constants were 1.864, 2.234, and 2.615 respectively for 
molecular nitrogen, molecular oxygen, and atomic oxygen. 


The neutral temperature was chosen as 1000 K at the base 
of the exosphere, the outer fringe region of the atmosphere. 
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The neutral temperature and densities as a function 
of altitude between 120 km and 250 km are shown in Figures 
3.4 and 3.5. There is still considerable uncertainty in 
determining the composition of the upper atmosphere. When 
comparing the model atmosphere adopted for this calculation 
and the model (S = 150) derived by Harris and Priester 
(1962), it can be seen that the density profile of each 
neutral constituent is relatively similar. The main difference 
between these two models arises from the different boundary 
conditions chosen for 120 km. These conditions are subject 
to change at different times of the year and at different 
places around the world. 

Energy transfer of electrons in elastic collisions 
with neutral particles in the upper atmosphere has been 
discussed by Hanson and Johnson (1961), Hanson (1963), and 
Dalgarno et al. (1963). The rate at which electrons lose 
their thermal energy to ions and neutral gases can be 
determined from the results derived by Desloge in 1962. 

The rate of exchange of kinetic energy between electrons 
and ions or neutrals can be expressed as following; 


dU 

dt 


m 8kT , 

-4 n n. - g - k ( — ^-) 2 Q_ ( T - t. ) 
elm. v m ' W D v e 1' 

1 e 


(3.7) 
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where n^, m^, and are the number density, mass and the 
Maxwellian temperature of the ion or neutral, k is the 
Boltzmann constant and Q D is the average momentum transfer 
cross-section, which can be expressed ass 


m 

_ e 

Q D = ( 2k T 


CO 

r 


m v 
e 




v q D (v) exp(- 


•) dv 


(3.8) 


where v is the electron velocity and q D (v) is the velocity 
dependent momentum transfer cross-section. Recently, Banks 
(1966) presented several energy loss rates for the individual 
gas components in the ionosphere derived from the available 
elastic cross-section data. His equation for the energy 
transfer rate due to molecular nitrogen was derived from 
the momentum transfer cross-section given by Englehardt 
et al. (1964), as follows* 

1.77 x 10 n e n(N 2 ) (1 - 1.21 x 10 T e ) T e ( T e ~ T n ) 
ev/cm 3 - 


(SO) 

^dt ; N, 


sec 


(3.9) 
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For the oxygen molecule, the adopted value of 
moment um transfer cross-section was based on Hake and Phelps's 
(1967) analysis. 


<S:> 0 2 = !- 21 x 10 ' 18 n e n ( 0 2 )(1 + 3 * 6 x 10 ‘ 2 T e )T e‘ T e ' V 

3 (3.10) 

ev/cm - sec 

For atomic oxygen, there are several experimental re- 
sults for the values of the total scattering cros-section, 
while none has been measured for the momentum transfer cross- 
section due to the chemical activity of oxygen atoms. Therefore 
the momentum transfer cross-section adopted for the atomic 
oxygen energy transfer rate was derived from the results of 
theoretical phase shift calculations which were done by 
Cooper and Martin (1962). The energy transfer rate can 
be expressed as; 


<f>0 = - 3 - 74 x 10 ' 18 n e nl0)T e (T e 


T ) ev/cm' 
n 


sec 


(3.11) 


In a recent paper, Dalgarno and Degges (1968) 
have shown that an efficient way of cooling the electron 
gas in the E and F regions of the ionosphere is by the 
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excitation of the fine structure levels of atomic oxygen 
through inelastic collisions with thermal electrons , as 
in the following equation: 

3 3 

e + 0( P T ) > e + 0( p_. ) 

J J 

the 3 P Q and 3 P^ levels lie 0.028 ev and 0.020 ev above 

3 

the P£ state. The importance of the heavy neutral 
particle analogue of equation (3.12) 


x + 0( 3 Pj) » x + 0( 3 P J( ) 


(3.13) 


in the thermal balance of the neutral atmosphere has been 
discussed by Bates (1951). 

The energy transfer rate shown below was calculated 
by Herman and Chandra (1969), using numerical values of 
the fine structure transition cross-section provided by 
Breig and Lin (1966). 


( dU) 

'“dt'O 


A 


- (3.7 - 2.92 x 10 4 T e ) ( 
x n e n(0) (T 0 - T n )/ T n 


9.06 + 6.57 x 10 4 T n ) 


x 10 


-13 


ev/cm^ - 


sec 


(3.14) 
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With increasing temperatures, cooling by 
excitation of the metastable 1 d level of atomic 
oxygen (Rees, et al., 1967) becomes significant. 
Secondary electron excitation predominates this process 
up to about 250 km. The rate of this energy transfer 
has been derived from the excitation cross-section of 
Smith et al., (1967), and is shown in the following 
equation : 


( §£) 0w = - 1.1 x 10 10 n e n(0)T^ e xp( 

(0.406 + 0.357 x 10 _4 T ) 

- (0.333 + 0.183 x 10' 4 T ) 

- (0.456 + 0.174 x 10" 4 T ) 

0 


- 2.27 x 10 4 /T ) 
0 


X exp(- 1.37 x 10 4 /T e ) 

X exp(- 2.97 x 10 4 /T ) 

© 


(3.15) 
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D. Results and Discussion 

The calculated rates of electron cooling by neutral 
constituents and ions as a function of altitude are tabu- 
lated in Table 3.3. The total electron cooling rate is 
the sum of electron cooling rates caused by ions and neutrals. 
The results are shown in Figure 3.6. The energy loss rates 
for six individual processes are shown in Figures 3.7, 3.8, 
3.9, 3.10, 3.11 and 3.12. From these tabulations, it is 
clear that the three dominant loss processes ares (1) ex- 
citation of fine structure levels of atomic oxygen, (2) 
elastic collision with and (3) electron-ion coulomb 

collision. The percentage losses for these three processes 
are plotted in Figure 3.13. 

A comparison of Figures 3.4 and 3.5 with Figure 3.6 
shows that the thermal electron cooling rate is approximately 
equal to the rate of heat input transfered from super- 
thermal electrons to thermal electrons. This fact bears 
out the assumption of local heating within the altitude 
ranges included in these calculations. 

Above 180 km, the cooling rate approximates the results 
of Butler and Buckingham (1962), which were derived from 
calculations based on test charged particles moving in a 
thermal plasma. Below 180 km, the cooling rate shows 
a closer agreement with Desloge's calculations (1962), 
which were derived by using two Maxwellian gases with 
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Fig. 3.9 THERMAL ELECTRON COOLING RATE BY 
ELASTIC COLLISION WITH 0. 



ALTITUDE (km) 



Fig. 3.10 THERMAL ELECTRON COOLING RATE BY 
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TABLE 3.3: ELECTRON COOLING RATES BY 

NEUTRAL CONSTITUENTS AND IONS 


Altitude 

Electron 

Cooling 

(km) 

f2 

°_2 

126. 

7,990 

763.0 

138. 

1,850 

147.0 

149. 

1,400 

99.1 

164. 

726 

44.4 

174. 

684 

38.3 

187. 

599 

32.1 

194. 

679 

35.0 

211. 

617 

30.7 

217. 

579 

28.5 

222. 

494 

23.4 

228. 

436 

20.2 

235. 

473 

21.1 

237. 

493 

21.9 

239. 

508 

21.4 

151. 

1,550 

105.0 

162. 

855 

52.5 

171. 

1,010 

60.5 

184. 

697 

37.7 

193. 

703 

36.9 

204. 

482 

23.9 

210. 

521 

25.7 

216. 

434 

21.0 

223. 

466 

22.0 

227 

393 

17.8 

231. 

. 419 

18.8 

235. 

439 

19.4 

237. 

463 

20.2 

239. 

513 

21.7 


Rates 

0 

(ev/cm^ - sec) 
0, 

O w 


A 

M 

494. 

22,000 

0.83 

165. 

6,230 

0.24 

150. 

4,490 

0.26 

103. 

2,700 

0.32 

114. 

2,660 

1.45 

129. 

2,740 

5.66 

163. 

3,140 

23.0 

209. 

3,600 

93.0 

218. 

3,640 

125.0 

207. 

3,530 

87.3 

201. 

3,480 

67.5 

239. 

3,940 

139.0 

253. 

4,140 

164.0 

265. 

4,280 

199.0 

169. 

4,800 

0.63 

120. 

3,150 

0.47 

164. 

3,640 

6.30 

144. 

3,020 

7.72 

175. 

3,370 

24.6 

143. 

2,760 

15.6 

172 

3,070 

47.9 

160 

2,930 

29.5 

194 

3,380 

63.6 

177 

3,170 

38.9 

197 

3,440 

58.9 

222 

3,760 

99.6 

242 

4,010 

132.0 

268 

4,310 

209.0 


+ 

0 


267. 

284. 

390. 

502. 

679. 

983. 

1280. 

2090. 

2410. 

3040. 

3340. 

4040. 

4450. 

4670. 

446. 

508. 

703. 

935. 

1260. 

1560. 

1770. 

2230. 

2930. 

3200. 

3700. 

4010. 

4430. 

4680. 
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different temperatures. These facts may suggest that 
the secondary electron distribution in the altitudes 
below 180 km are closer to a Maxwellian distribution 
than distributions of the secondary electrons above 
180 km. 

In calculating electron energy losses to neutral 
constituents, it is necessary to take into account 
the uncertainties existing in the neutral density 
measurements in the atmosphere. Since neutral densities 
vary at different places around the world and at different 
times of the year, it is difficult to choose a model neutral 
densities profile which represents the atmospheric density 
prevailing at Wallops Island, Virginia, when the exper- 
iment was performed. Recently Hedin et al. (1970) have 
pointed out the fluctuations up to 20 percent of neutral 
densities can exist in the atmosphere and that fluctuations 
of atmospheric density may be related to the gravity waves 
in the upper atmosphere. Gravity waves are very low 
frequency atmospheric waves. Mathematically, the 
existance of these waves is derived from gravity terms 
in the equation of motion. At present, the physical 
sources of these waves are still not known. They may be: 

(1) turbulence in the lower atmosphere, or (2) pertur- 
bations of the wind caused by mountains. Other possible 
sources may be thunderstorm activity and nuclear 
explosions. 
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In this investigation, the calculated electron energy 
loss rate was largely affected by the losses due to exci- 
tation of the fine structure levels of atomic oxygen. 
Therefore, the number density of atomic oxygen chosen for 
calculating the energy loss rate would have affected 
significantly the accuracy of the calculated energy balance. 
Due to the fluctuation of neutral density as pointed out 
by Hedin et al. (1970), density fluctuation could have 
caused an error of approximately 10 percent in the calculated 
energy loss rate. The energy losses due to elastic collision 
of electrons with N 2 and 0 2 represented only a small portion 
of the total energy losses; therefore, the variation of N 2 
and C >2 densities in the upper atmosphere could have affected 
the calculated energy loss rate only slightly. 

The energy loss by vibrational excitation of nitrogen 
molecules has been studied by Dalgarno et al. (1968) 
following the. procedures described by Rees et al. (1967). 
According to their calculations, this process may be a 
dominant energy loss mechanism for electrons in the energy 
interval of roughly 1.7 to 3.5 ev; therefore, the energy 
distribution function in this region should show a strong 
perturbation. Hoegy et al. (1965) and Shea et al. (1968) 
included this vibrational loss process in deriving the 
energy distribution function for photoelectrons and showed 
a deep valley in the distribution function around 2 ev 
region. Since their calculations were based on the momentum 
transfer cross-section measured in the laboratory, their 
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opinions on vibrational energy loss in the upper atmosphere 
are a matter of controversy. In this experiment, the retarded 
probe current collected showed two distinct linear regions 
(Figure 2.4) when the current was plotted on a logarithmic 
scale. If the energy losses due to vibrational excitation 
of molecular nitrogen were taken into account in this 
investigation, the total energy loss rate would increase 
by 100 percent as shown in Figure 3.14, and its effect 
upon the distribution function could be shown on the 
current collected by the probe. The data obtained by this 
experiment, however, did not show such perturbation in 
the energy distribution function between 1.7 and 3.5 ev as 
calculated by Hoegy et al. (1965) and Shea et al. (1968). 

As a result, it may be suggested that the vibrational energy 
loss process does not have a large effect on the electron 
energy distribution function in the altitude ranges between 
120 and 240 km. For this reason, the vibrational energy 
loss by molecular nitrogen has not been adopted in this 
dissertation. 

Estimated errors for energy losses vary according 
to the processes involved. The most important energy loss 
mechanism between 120 and 230 km is the excitation of fine 
structure levels of atomic oxygen, and the accuracy of the 
loss rate depends upon the cross-section adopted. Herman 
and Chandra (1970) have estimated that their calculated 
loss rate may vary approximately 5 percent. As a result, 



ALTITUDE (km) 



Fig. 3.14 THERMAL ELECTRON COOLING RATE BY 

EXCITATION OF VIBRATIONAL STATE OF N 2 
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an error of approximately 3 percent could occur in the 
calculated loss rates. Elastic collision of electron 
with molecular nitrogen at lower altitudes (~130 km) can 
cause an energy loss rate of approximately 30 percent. 

The uncertainty of the cross-section in this process is 
believed to be within 15 percent. An error of approximately 
5 percent could apply in this case. Electron-ion coulomb 
collision at higher altitudes (^235 km) dominate 50 percent 
of the total loss rate, the accuracy of the electron-ion 
coulomb collision cross-section is estimated to be 
approximately 10 percent (Banks, 1966). Therefore, the 
optimum estimate of error for the calculated total heat 
loss rate is 10 - 20 percent. Estimations of errors for 
•the individual processes have been calculated for the different 
altitudes, and the results are plotted in Figure 3.13. 
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CHAPTER IV 

The Production of Super thermal Electrons 
A. Introduction 

The primary objective of this investigation is to 
study the energy distributions of electrons in the energy 
range from 0 to 15 ev in the middle ionosphere. To 
accomplish this objective, the author developed a theoretical 
model of the electron energy distributions in the above 
energy range. In developing his model, the author did not 
accept the premise that the vibrational excitation of 
molecular nitrogen in the ionosphere represents an important 
mechanism of energy losses for electrons. In this study, 
superthermal electrons refer to those electrons with energies 
greater than ^ 1 ev, without particular reference to their 
source of origin. It is believed that the main sources 
of these superthermal electrons are primary photoelectrons 
and secondary electrons produced by the primary photoelectrons. 
At least some superthermal electrons are produced through 
superelastic collisions in which the energy of electrons 
is increased by interactions with either ions or neutrals 
in excited or metastable states. 

If it is assumed that photoelectrons are the dominant 
portion of superthermal electrons, then it can be expected 
that at high energies the superthermal electron energy 
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distributions would approximate the primary photoelectron 
distribution. Calculations for the energy distribution of 
photoelectrons have been done by Hoegy, et al. (1965), 

Shea, et al. (1968), and Dalgarno, et al. (1969). They 
all use the solar flux data provided by Hinteregger, et 
al. (1965) as the source for the production of photoelectrons, 
but have adopted different cross-sections and weighting 
factors for loss processes. The process of the production 
of secondary electrons has remained a problem. Green and 
Barth (1967) mentioned it in their calculations, but the 
problem has not been completely investigated. 

The equilibrium energy distribution function is expressed 

as i 


F (E) 


j p ( E 1 ) dE 1 

P ( E) _ _ E 

dE 

dt dt 


(4.1) 


where — - is the energy loss rate of the photoelectrons, 
and P(E) is the accumulated production rate defined by 


P ( E) 


| p(E« ) dE 1 


(4.2) 


E 
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and where p(E' ) dE 1 is the rate at which photoelectrons are 
produced in the energy interval E 1 and E 1 + dE' . In this 
chapter the equilibrium energy distribution will be cal- 
culated by using the above equation. Discussions will be 

made concerning the accumulated production rate P(E) and 

c3E 

the energy loss rate — . Also included will be a 
possible explaination of the linear section of the probe 
current described in Figure 2.4. In addition, some of the 
calculated results and discussions will be presented. 

B. Production of Photoelectrons 

Energetic electrons are produced in the upper atmosphere 
by absorption of extreme ultraviolet solar radiation (EUV). 

In thermospheric physics, the EUV refers mainly to the 

o o 

range of wavelengths from about 1750 A to 170 A, together 
with X-rays, including both soft and hard, but generally 
not referring to Y-rays. Photons with energies greater 
than about 12 ev can ionize one or more of the major 
atmospheric constituents are tabulated in Table 4.1. The 
energy of ejected photoelectrons depends on the energy 
of the incident photon and the ionization potential of the 
target atom. Tohmatsu, et al. estimated that about two- 
thirds of EUV energies are used in removing electrons 
from atmospheric molecules. The remaining one-third are 
transferred to the kinetic energy of the ejected photo- 


electrons 



TABLE 4.1 


Ionization Potentials for Photoionization of Atoms and Molecules 


First Ionization Potential (ev) 


NO 

9.25 

°2 

12.08 

H 

13.59 

0 

13.61 

N 

14.54 


15.58 

He 

24.56 
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In order to calculate the accumulated production rate 
P ( E ) for the primary photoelectrons, it is necessary to 
know the probability of ionization in each of the electronic 
states of ions in addition to the ionization cross-sections. 
Tohmatsu, et al. (1965) performed a calculation of the 
kinetic energy spectrum of the primary photoelectrons. 

The solar flux data, ionization cross-sections, and absorption 
cross-sections were compiled by Hinteregger, et al. (1965). 
Tohmatsu, et al. (1965) also calculated the photoelectron 
production rates for different solar zenith angles and for 
several different altitudes ranging from 150 km to 500 km. 

o 

For calculations in this study, the solar zenith angle ~ 30 
was used, and the accumulated photoelectron production rate 
can be written as: 


P(E) = P exp(-E/E ) (4.3) 

where P Q represents the total production rate, and E q is 
the mean kinetic energy of the photoelectrons. The mean 
energy of the primary photoelectrons is approximately 
constant at ^ 10 ev above 200 km. Below 200 km, the mean 
energy becomes large due to the increasing importance of 

O 

He II 304 A emission. The mean kinetic energy for the 
photoelectrons used in this calculation was taken from 
Tohmatsu, et al. (1965) and tabulated in Table 4.2. 



TABLE 4.2 


Total Production Rates and Average Kinetic Energy for Photoelectrons 


Altitude 

Total Production Rate 

Average Kinetic Energy 

(km) 

, -3, 

(cm ) 

(ev) 

150 

6,000 

15 

200 

2,000 

10 

240 

700 

10 
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Photoelectron production depends upon the EUV solar 
photons in the upper atmosphere. The total production 
rates as a function of altitude have been derived by 
Hinteregger, et al. (1965) and are reproduced in Figure 

4.1. The figures for the total production rate and mean 
kinetic energy for the photoelectrons used in the calculations 
in this study are shown in Table 4.2. By utilizing equation 
(4.3), the accumulated production -rate for the photoelectrons 
as a function of energy can be derived. The production 
rates for several different altitudes are plotted in Figure 

4.2. 

C. Cooling Processes of Photoelectrons 

Photoelectrons lose energy by exciting and ionizing 
the neutral particle constituents of the atmosphere and 
by elastic collisions with the ambient thermal electrons. 
Elastic collisions with neutral particles absorb a 
negligible ’portion of the photoelectron energy. 

Energy loss due to collision with ambient electrons 
can be readily estimated by equation (3.1). When the energies 
of the photoelectrons are substantially greater than those 
of thermal electrons, equation (3.1) can be simplified 
and written ass 

-4 _p 

= 1.16 x 10 n e E 2 ev/sec 


dE 

dt 


(4.4) 
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ACCUMULATED PRODUCTION RATE (cm -3 sec 
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where E is the fast electron energy in ev. 

The rate of energy loss of photoelectrons due to 
inelastic collisions with a neutral species n(m) is shown 
in the equation: 

= n(m) E Q d v ev/sec (4.5) 

where n(m) is the number density of the neutral species, 

E and v are the energy and the magnitude of the velocity 
of the photoelectron, and Q D is the momentum transfer 
cross-section. 

In measurements used in this study, the secondary 
energy distribution were devoid of the structure suggested 
by Hoegy, et al. (1965), and Shea, et al. (1968). The two 
groups suggested a depression in the distribution curves 
at about 3 ev is related to the vibrational excitations of 
N 2 , and a simultaneous decrease in vibrational and electronic 
excitation energy losses should be the cause of the hump 
between 5 and 6 ev. Due to the uncertainties in determining 
the inelastic collision cross-sections for N 2 , C> 2 , and 0, 
it is extremely difficult to adopt a set of reliable cross- 
sections. Calculations used in this study depict a smooth 
photoelectron energy distribution curve devoid of structure 
based on the energy loss rate of electrons to the neutrals 
as presented by Shea, et al. (1968). 
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The calculated energy loss rates of superthermal electrons 
due to inelastic collisions with neutral particles and 
elastic collisions with thermal electrons are shown in 
Figures 4.3 and 4.4. The total time rates of energy loss 
are given in Figure 4.5. From these curves it can be 
seen that at an altitude of 150 km the energy loss to the 
neutrals is more important than the electron-electron loss. 

At 240 km the electron-electron energy loss becomes the 
dominant loss for electron energy below 5 ev. Above 5 
ev, the inelastic losses always are the important loss 
for photoelectrons. 

D. Equilibrium Energy Distributions of Superthermal Electrons 
Apart from the thermal electron gas, electrons with 
energies greater than approximately 1 ev are called 
"superthermal electrons." In this study it was assumed 
that they were originally all primary photoelectrons. 

Based on this assumption, a formal expression for the equili- 
brium energy distribution function can be constructed as 
follow: 


p(E)dE = F(E)|| / E - F(E)§| / E+dE = 


- d 


F ( E) 


dE 

dt 


(4.6) 


or 
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oo 


F ( E) . = 


p(E* ) dE ' 


dE 

dt 


(4.7) 


Both the accumulated production rate and the energy 
loss rate have been calculated in Sections B and C of this 
chapter. Hence, by using equation (4.7) the distribution 
function F(E) can be obtained immediately. It has been 
calculated for several different altitudes, and the results 
are shown in Figure 4.6. 

Comparisons of the electron energy distribution 
function published by Shea, et al. for 156 km (1968), 

Hoegy, et al. and Nagy and Fournier for 160 km (both 1965), 
and the calculated results for 150 km are plotted in 
Figure 4.7. The main difference between the calculated 
results achieved in this study and the results published 
by the two groups above are indicated by the existence of 
the trough around 3 ev and the hump around 5 ev, due to 
vibrational loss processes, of the energy distribution 
function. 

In Section A of Chapter II, it was shown that the 
energy distribution function F(E) can also be evaluate by 
taking the second derivative of the probe current 
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(Druyvesteyn relationship). It is difficult to maintain 
accuracy in performing the graphical differentiation without 
smoothing the raw data. Therefore no attempt was made to 
obtain the distribution function F(E) by using this method. 

E. Delta Function in the Distributions 

In a typical section of the i-V curve, as shown in 
Figure 2.5, the current appears to be relatively linear 
with applied retarding potential in the range of 6 - 13.5 
volts. Since the second derivative of the current with 
respect to the applied voltage is nearly zero in this region, 
according to the Druyvesteyn relationship given in Section A 
of Chapter II, it may be concluded that relatively few 
electrons exist in the energy interval between 8 and 15 ev. 

The recorded current in this range of retarding potential 
could be due mainly to electrons with energies greater 
than 15 ev. In order to produce a linear current in the 
energy interval between 8 and 15 ev, a particular distribution 
function, which is different from the exponential form is 
needed. It can be assumed that electrons with energies 
greater than ~ 15 ev are approximately monoenergetic with 
energy E q . Based on this assumption, the current due to 
the delta function can be derived. 

By using the notations defined in Section B of Chapter II, 
the current related to the delta function may be written as: 
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i 1 = 2 ^ale J dv / v^f(v^,v) dv^ 

“ Vma/ oi 

V max . 7(f) 2 (V. 2 - ■ % *■ «- <* 

velocity of the monoenergetic particles with energy E„ , 
and the distribution function can be expressed as: 


/ 2 2 

r 2 / , V P + V x 2 2 / 2 

■ n 0 /(27 Tv 0 /I 5 — ) , v + v S v o 

f (v v) = ) J V 0 2 r 

[0 v r + v > v 


and satisfied the normalization condition 


Jc 2 - v; 


dv I dv, f(v^ ,v) = n 0 


"Jc 2 - ^ 


and where n c is the particle density for electron with 
energy E c . 

The current i 1 due to the delta distribution can be 
derived by substituting equation (4.9) into equation (4.8). 


Vinax /- J c z - V 2 I 

= 27raln c e dv I v p /(27 Tv 0 2 /l - 


v p 2 ♦ V 2 

— > av f 

v o r 


irRln 0 ev 0 


11 - lr ) 


(4.12) 
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The current i' due to the delta function in the equilibrium 

distribution is a function of the probe radius and the 

applied potential, and also depends on the value of E q . 

For a given E q , it is possible to calculate the current 

as a function of the applied retarding potential for a 

given probe with radius R. The relation of the current and 

the retarding probe potential has been calculated for probe 

radius equals to 1 cm, with E q equals to 25 ev, and the 

results are shown in E'igure 4.8. 

F. Results and Discussion 

From Figure 4.6 it can be seen that above 3 ev the 

distribution function F(E) is an exponentially decreasing 

function of energy. By integrating the area below the 

distribution curves, it should be possible to calculate 

approximately the superthermai electron density. At 150 km 

and 200 km, the integrated values are approximately 100 

electrons/cm^. These values are comparable to the measured 

superthermai electron concentrations. At 240 km, the 

3 

calculated value is approximately ~ 50 electrons/cm , while 

the measured density is 4 to 5 times larger. This fact may 

suggest that photoionization by incident solar radiation 

is not the only source for these energetic electrons. 

The concentration of electrons with energies greater 

than 15 ev increases from ^ 5 elect rons/cnf 5 at altitudes 

3 

of approximately 150 km to — ' 15 electrons/cm at about 
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240 km. Figure 4.6 shows that the energy distribution 
function for 150 km is larger than that for 240 km. An 
extrapolation would indicate that the concentration of 
electrons with energies greater than 15 ev would not be 
larger at 240 km than at 150 km. This difficulty can be 
resolved by superimposing the delta function on the energy 
distribution. By choosing different electron densities for 
electrons with 25 ev energy at 150 km and at 240 km, this 
problem can be avoided. 

The uncertainty in the measurements of the EUV solar 
spectrum will also affect the calculations significantly. 

The solar spectrum measurements used in this study which 
were taken by Hinteregger, et al. in (1965), were subjected 
to variations in place and time. Several Orbital Solar 
Observatory satellites put into orbit since 1962 should 
eventually provide better data on the solar EUV spectrum 
measurements,' although, up to the present time no data have 
been released. 

It has been shown that when the ejected primary photo- 
electron has a kinetic energy greater than the ionization 
potentials of the atoms or molecules it encounters, the 
production of secondary electrons occurs. When electrons 
have energies greater than 70 ev, they produce one ion-pair 
for every ^ 34 ev. Electrons with energies between 70 ev 
and 20 ev may produce one ion-pair. Dalgarno, et al. (1963) 
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indicated that at approximately 120 km the majority of the 
photoelectrons being produced have energies between 60 and 
120 ev. Above 150 km, the production spectrum is dominated 
by photoelectrons produced in the energy interval of 25 to 
35 ev. Therefore, it can be assumed that the secondary 
electrons produced by primary photoelectrons may contribute 
substantially to the total electron production. Theoretical 
considerations concerning the energy distributions of the 
secondary electrons must begin with a careful examination 
of the relevant cross-sections of the production processes. 

This examination should include the study of impact ionizations, 
superelastic collisions, and inelastic collisions. Experimental 
data are available concerning impact ionization and inelastic 
collisions by electrons for N 2 , 0 2> and O, but these 
measured cross-sections usually are not very accurate. 

Under these conditions, the secondary electron distribution 
requires further investigation. 

Figure 4.7 shows that a monoenergetic distribution will 
produce a linear section of current on the i-V plot. In the 
altitude ranges used in this investigation, the particle- 
particle interactions do not seem to be the cause of the 
delta function in the energy distributions. One possible 
reason for the existence of a monoenergetic distribution is 
the interaction of wave-particle in the ionosphere (Sinha, 

1970). The emerging energetic electrons interact with electro- 
magnetic waves in the ionosphere to form a steady peak in 
the energy distribution. 



86 


CHAPTER V 

Summary and Conclusions 


A. Summary 

1. Results of Calculations 

The energy balance of electrons in the middle 

ionosphere can be calculated by using measured superthermal 

and thermal electron densities and temperatures, in combination 

with the model atmosphere and electron collision cross- 

sections. In the calculations used in this study, the heat 

input rates for superthermal electrons to thermal electrons 

3 3 

are found to be in the range between 3 x 10 ev/cm -sec and 
4 3 

1.5 x 10 ev/cm -sec. The numbers in this range correlate 
closely with the calculated total heat loss rate for electrons 
to ions and neutral constituents. The deviation for heat 
loss and input becomes noticeable to approximately 235 km, 
which is close to the apogee (239 km) for the rocket flight 
which provided data used in this investigation. 

In the rocket flight, the measured superthermal 
electron densities were approximately 150 - 200 electrons/ 

3 

cm of heights between 150 and 200 km. In this altitude 
range, the number density of superthermal electrons corresponded 
to 0.02 percent of the thermal electron density. This 
may suggest that superthermal electrons at varying heights 
in this altitude range were produced by similar physical 
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processes. By assuming a photoionization origin for 
supertbermal electrons, it was possible to calculate the 
equilibrium energy distribution function. 

For electrons with thermal energy, electron-ion 
collisions and elastic collisions with neutral constituents 
are the dominant energy loss processes. Below 300 km, 
the elastic collisions constitute the dominant energy loss 
processes for thermal electrons. Inelastic collisions and 
ionization are the most important energy loss processes 
for the superthermal electrons throughout the altitude range 
of this study. 

2. Possible Errors in the Calculations 

The overall correlation between the heat input 
rate for superthermal electrons to thermal electrons and 
the heat loss rate for thermal electrons to neutral constituents 
and ions is gratifying. A large number of theoretical studies 
and experimental works have been carried out to improve 
the accuracy of collision cross-sections. Most of the 
collision cross-sections concerning processes in the 
ionosphere are measured under laboratory conditions which 
are far different from those found in the ionosphere. For 
example, the measurements of the momentum transfer cross- 
section for molecular nitrogen were performed in a container 
filled with pure nitrogen gases, while in the ionosphere 
besides nitrogen molecules there exist oxygen atoms and 
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oxygen molecules. Recently, Breshears and Bird (1968) 
reported in an experiment which measured the vibrational 
relaxation time for molecular nitrogen that the existence 
of atomic oxygens would decrease the relaxation time by 
nearly two order of magnitude, in a comparison of values 
measured in pure nitrogen. This phenomenon shows that the 
physical parameters in the ionosphere may be considerably 
different from those measured in laboratory experiments. 

For a better understanding of the energy balance in the 
ionosphere, it is desirable to measure collision cross- 
sections under ionospheric conditions. 

Another source of possible experimental error 
may arise if a component of the measured flux consisted 
of photoelectrons originating on surfaces of rockets. It 
is extremely difficult to evaluate these photoelectrons 
quantitatively. In Huang's experiment, the pulse probe 
was mounted ~ 20 in. above the main payload section. Hence, 
the rocket electrons would require a large angle scatter 
to be collected by the pulse probe. The probability that 
electrons undergo a large angle scattering should be small, 
and can be neglected in this calculation. The average 
energy for electrons from an aluminum surface exposed to 
sunlight in the upper atmosphere is approximately 0.2 ev, 
with electrons seeming not to affect the concentration of 
measured superthermal electrons. In the experiment, the 
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photoelectrons from the collector were suppressed by a 
large positive potential bias on the collector and the 
innermost grid; therefore, they should not have produced 
an appreciable photoelectron current. 

Photoelectrons with several ev energies can travel 
a large distance before they finally become thermalized. 

Recently, Heikkila and Winningham (1970) reported that a 

8 2 

flux of 10 electrons/cm -sec for escaping electrons was 
measured in the midlatitude daytime ionosphere. Their 
measurements confirmed the fact that the fast electrons 
could travel a large distance and did not deposit their 
energy locally. Therefore, this nonlocal heating effect 
may have had some influence on the results of this study. 

No quantitative calculations concerning the importance 
of this effect in the altitude range between 120 and 240 km 
have been performed. Further study of this problem is 
highly recommended. 

B. Conclusions 

The successful pulse probe experiment used in 
this study detected both the thermal and the secondary 
electron distributions. These measured values provided 
an opportunity to calculate the equilibrium distribution 
of electrons for energies up to 15 ev. To determine the 
electron energy distribution beyond 15 ev, the applied 
retarding potential for the probe should be increased to 
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about 50 ev. Under this circumstance, the calculation of 
the equilibrium distribution of electrons up to approximately 
50 ev can be performed, giving a better understanding of 
the secondary electron distributions. 

The calculated results from measurements taken 
during the successful pulse probe experiment indicate 
the followings 

1) The total heat input rate of superthermal 

3 3 

to thermal electrons fluctuates between 3 x 10 ev/cm -sec 

4 3 

and 1.5 x 10 ev/cm -sec in varying altitudes between 120 
and 240 km. These values correspond to the calculated 
thermal electron heat loss rate. 

2) The electron-ion Coulomb collisions dominate 
the energy loss for thermal electrons above 200 km. The 
excitation of fine structure levels of atomic oxygen is 
also an important energy loss process in the altitude 
range of thi-s study. 

3) The equilibrium distributions for electrons 
constitute an exponential function for electron energies 

up to 15 ev. 

4) The linear probe current corresponding to 
the retarding potential in the interval 8 to 15 volts 
can be reproduced by a monoenergetic energy distribution 


at '■V' 25 ev 



91 


A satisfactory description of equilibrium distributions 
for electrons with energies below 15 ev has been derived 
through a combination of information concerning atmospheric 
compositions, solar flux, and cross-sectional data. The 
theoretical model used in this study assumed that production 
and energy losses of photoelectrons occurred locally. More 
sophisticated calculations concerning production and energy 
losses of photoelectrons in the middle ionosphere could be 
made possible through the development of a theoretical model 
which would describe nonlocal electron heating processes. 
Hopefully this type of model would provide desirable data 
concerning electron energy balances at altitudes below 250 km. 
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